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ROAD SURFACE FRICTION STATE ESTIMATING APPARATUS AND 

METHOD THEROF 

INCORPORATION BY REFERENCE 

[0001] The disclosure of Japanese Patent Application No. 2002-234505 
filed on August 12, 2002 including the specification, drawings and abstract is 
incorporated herein by reference in its entirety. 

BACKGROUND OF THE INVENTION 

1 . Field of the Invention 

[0002] The invention relates to a road surface friction state estimating 
apparatus. More particularly, the invention relates to a road surface friction state 
estimating apparatus which estimates a road surface friction state using self aligning 
torque generated in a tire of a vehicle, and a method thereof. 

2. Description of the Related Art 

[0003] As related art 1 of the invention, Japanese Patent Laid-Open 
Publication No. 6-221968 discloses a road surface friction coefficient detecting 
apparatus which estimates a grip state using self aligning torque (hereinafter, referred 
to as "SAT"). 

[0004] The apparatus according to the related art 1 estimates a grip state 
based on SAT and a front wheel lateral force which are calculated based on vehicle 
movement. According to the related art 1, it is possible to accurately estimate a grip 
state even when the vehicle runs on a bank road surface, or when a relation between a 
lateral force and a slip angle is non-linear. 

[0005] However, the apparatus according to the related art 1 is likely to be 
influenced by road surface disturbance. Also, when steering is performed quickly, the 
phase of the calculated SAT value and that of the calculated lateral force value are 
deviated from each other. Accordingly, a grip state cannot be accurately estimated. 

[0006] As another related art 2 of the invention, Japanese Patent 
Application No. 2001-212683 discloses art in which a front wheel slip angle is 
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estimated based on a steering angle, SAT is estimated by removing friction of a 
steering system from a sum of steering torque supplied by a driver and assist torque 
applied to a power steering device, and a grip state is estimated based on the front 
wheel slip angle and SAT. 
5 [0007] However, according to the related art 2, a grip state cannot be 

accurately estimated when the relation between the lateral force and the slip angle 
becomes non-linear and accordingly the front wheel slip angle cannot be accurately 
estimated in the case where the steering angle is 0 and SAT is not 0 during running on 
a bank road surface. 

10 

SUMMARY OF THE INVENTION 



[0008] In order to solve the aforementioned problem, the invention is 
proposed. Accordingly, it is an object of the invention to provide a road surface 

15 friction state estimating apparatus which can accurately estimate a road surface 

friction state without being influenced by road surface disturbance even when steering 
is performed quickly, and a method thereof. 

[0009] A road surface friction state estimating apparatus according to a 
first aspect of the invention includes a self aligning torque estimating portion which 

20 estimates self aligning torque; a slip angle estimating portion which estimates a slip 
angle; a high pass filter which performs high pass processing on the slip angle 
estimated by the slip angle estimating portion; a lateral force calculating portion 
which calculates a lateral force; a slip angle converting portion which converts the 
lateral force calculated by the lateral force calculating portion into a slip angle; a low 

25 pass filter which performs low pass filter processing on the slip angle obtained by 
conversion by the slip angle converting portion; a summing portion which sums the 
slip angle on which the high pass filter processing has been performed by the high 
pass filter and the slip angle on which the low pass filter processing has been 
performed by the low pass filter; and a road surface friction state estimating portion 

30 which estimates a road surface friction state based on the slip angle obtained by 

summation by the summing portion and the self aligning torque estimated by the self 
aligning torque portion. 

[0010] The self aligning torque estimating portion estimates self aligning 
torque generated in a tire. The invention is not limited to a specific method of 
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estimating self aligning torque. The slip angle estimating portion estimates the slip 
angle of the tire. 

[0011] The high pass filter performs high pass filter processing on the slip 
angle estimated by the slip angle estimating portion, thereby removing a drift error 
5 contained in the slip angle estimated by the slip angle estimating portion, and 

extracting a high-frequency component of the slip angle, which does not have a phase 
lag with respect to the self aligning torque. 

[0012] The lateral force calculating portion calculates the lateral force 
generated in the tire. Since there is a substantially linear relation between the lateral 

10 force and the slip angle, the slip angle converting portion converts the lateral force 

into the slip angle based on this relation. The low pass filter performs low pass filter 
processing on the slip angle obtained by conversion by the slip angle converting 
portion, thereby removing a variation component such as disturbance noise contained 
in a high frequency region, and extracting an accurate low-frequency component 

15 during running on a bank road. 

[0013] The summing portion sums the slip angle on which the high pass 
filter processing has been performed and the slip angle on which the low pass filter 
processing has been performed, thereby obtaining the slip angle which does not have 
disturbance noise, a drift error, or a phase lag with respect to the self aligning torque. 

20 [0014] The road surface friction state estimating portion estimates a road 

surface friction state, that is, a friction state between the tire and the road surface, 
based on the slip angle obtained by summation by the summing portion and the self 
aligning torque estimated by the self aligning torque estimating portion. 

[0015] Therefore, according to the first aspect of the invention, it is 

25 possible to accurately estimate a road surface friction state without being influenced 
by disturbance noise or a drift error even during running on a straight horizontal road 
or a bank road. 

[0016] Further, in the road surface friction state estimating apparatus 
according to the first aspect of the invention, the high pass filter and the low pass 
30 filter may have the same breakpoint frequency. Further, the sum of a transfer function 
of the high pass filter and a transfer function of the low pass filter is preferably 1. 
Thus, when the same signal is input to the high pass filter and the low pass filter, and 
outputs of both the filters are summed, the original signal is restored. 
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[0017] Therefore, according to the aspect, it is possible to obtain a slip 
angle which does not have disturbance noise or a drift error. Thus, it is possible to 
accurately estimate a road surface friction state using the slip angle. 

[0018] Also, the value of the breakpoint frequency of the high pass filter 
5 and the low pass filter may be set according to a vehicle movement state. The vehicle 
movement state includes a vehicle speed, a turning state, and the like. 

[0019] For example, when a vehicle enters a bank road from a straight 
horizontal road, the breakpoint frequency of the high pass filter and the low pass filter 
needs to be equal to or higher than a road surface cant change speed, that is, a speed at 
10 which a cant of a road surface changes. The road surface cant change speed is 

proportional to the vehicle speed. Therefore, when the breakpoint frequency is set so 
as to be higher in proportion to the vehicle speed, it is possible to accurately estimate 
a road surface friction state even when the vehicle enters a bank road. 

[0020] For example, when a spin occurs, the breakpoint frequency of the 
15 high pass filter and the low pass filter may be set so as to be higher as a deviation 
between the slip angle estimated by the slip angle estimating portion and the slip 
angle obtained by conversion by the slip angle converting portion is larger. 

[0021] The deviation is large when the vehicle runs on a bank road, or 
when the relation between the lateral force and the slip angle is non-linear. In such a 
20 case, the accuracy of the slip angle estimated by the slip angle estimating portion 
deteriorates. Therefore, it is preferable to use the slip angle into which the lateral 
force is converted. Thus, when the breakpoint frequency is set in the manner 
described above, it is possible to accurately estimate a road surface friction state even 
if the vehicle suddenly enters a bank road, or a spin suddenly occurs. 
25 [0022] Accordingly, by setting the breakpoint frequency according to a 

vehicle movement state, it is possible to accurately estimate a road surface friction 
state even if the vehicle enters a bank road from a straight horizontal road, or a spin 
suddenly occurs. 

[0023] Also, in the road surface friction state estimating apparatus 
30 according to the first aspect of the invention, the road surface friction state estimating 
portion may calculate a self aligning torque reference value based on the slip angle 
obtained by summation by the summing portion, and may estimate a grip degree 
based on a ratio between the self aligning torque reference value and the self aligning 
torque estimated by the self aligning torque estimating portion. 
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[0024] Thus, according to the aspect of the invention, the self aligning 
torque reference value is calculated based on the slip angle which is not influenced by 
road surface disturbance and which does not have a drift error, and the grip degree is 
estimated based on the ratio between the self aligning torque reference value and the 
5 self aligning torque. Therefore, it is possible to accurately estimate the grip degree 
equivalent to an excess of a frictional force in the lateral direction. 

[0025] Further, the road surface friction state estimating portion may 
estimate a road surface friction coefficient based on the grip degree and the lateral 
acceleration. 

10 [0026] By using the estimated grip degree when estimating the road 

surface friction coefficient, it is possible to accurately estimate the road surface 
friction coefficient. The road surface friction coefficient may be estimated when the 
grip degree is equal to or lower than a predetermined reference value. 

[0027] A method of estimating a road surface friction state according to a 

15 second aspect of the invention includes the steps of estimating self aligning torque; 
estimating a slip angle; performing high pass filter processing on the estimated slip 
angle; calculating a lateral force applied to a tire; converting the calculated lateral 
force into a slip angle; performing low pass filter processing on the slip angle 
obtained by conversion; summing the slip angle on which the high pass filter 

20 processing has been performed and the slip angle on which the low pass filter has 
been performed; and estimating a road surface friction state based on the slip angle 
obtained by summation and the estimated self aligning torque. 

[0028] According to the second aspect of the invention, it is possible to 
obtain the slip angle which does not have disturbance noise or a drift error. Thus, by 

25 using the slip angle, it is possible to accurately estimate a road surface friction state. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0029] The foregoing and further objects, features and advantages of the 
30 invention will become apparent from the following description of preferred 

embodiments with reference to the accompanying drawings, wherein like numerals 
are used to represent like elements and wherein: 
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FIG. 1 is a block diagram showing a configuration of a road 
surface friction state estimating apparatus according to an embodiment of the 
invention; 

FIG. 2 is a block diagram showing a functional configuration of an 
5 electric control unit (ECU); 

FIG. 3 is a chart showing an estimated SAT value with respect to a 
sum of steering torque and assist torque; 

FIG. 4 is a chart showing the estimated SAT value with respect to 
the sum of the steering torque and the assist torque, for describing a method of 
10 removing a hysteresis characteristic; 

FIG. 5A is a chart showing a time-dependent change in the sum of 
the steering torque and the assist torque during running on a road having a high |i 
value, and FIG. 5B is a chart showing a time-dependent change in the estimated SAT 
value during running on the road having a high ja value; 
15 FIG. 6 A is a chart showing a relation of the sum of the steering 

torque and the assist torque with respect to a slip angle during running on a road 
having a high \i value and during running on a road having a low ja value; and FIG. 
6B is a chart showing a relation of the estimated SAT value with respect to the slip 
angle during running on the road having a high [i value and during running on the 
20 road having a low |i value; 

FIG. 7 is a chart showing time-dependent waveforms of an 
estimated slip angle value a E , a converted slip angle value ot T , and an integrated slip 
angle a! when the vehicle enters a bank road; 

FIG. 8 is a chart showing the estimated SAT value with respect to 
25 the slip angle during running on a dry asphalt bank road; 

FIG. 9 A is a chart showing a grip degree during running on a dry 
asphalt bank road, and FIG. 9B is a chart showing a time-dependent waveform of the 
steering angle; 

FIG. 1 OA is a chart showing the grip degree during running on an 
30 ice road at a speed of 40 km/h, FIG. 10B is a chart showing a road surface |i est value 
(an estimated \i value), and FIG. 10C is a chart showing a time-dependent waveform 
of the steering angle; 

FIG. 1 1 A is a chart showing the grip degree during running on a 
pressed snow road at a speed of 40 km/h, FIG. 1 IB is a chart showing the road 



surface ^ est value (the estimated |x value), and FIG. 1 1C is a chart showing a time- 
dependent waveform of the steering angle; and 

FIG. 12A is a chart showing the grip degree during running on an 
asphalt road at a speed of 40 km/h, FIG. 12B is a chart showing the road surface ja est 
5 value (the estimated p, value), and FIG. 12C is a chart showing a time-dependent 
waveform of the steering angle. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

10 [0030] Hereinafter, preferred embodiments of the invention will be 

described in detail with reference to the accompanying drawings. 

[0031] FIG. 1 is a block diagram showing a configuration of a road 
surface friction state estimating apparatus according to an embodiment of the 
invention. The road surface friction state estimating apparatus can be employed, for 

15 example, in a vehicle provided with an electric power steering device. Also, the road 
surface friction state estimating apparatus can be employed in a vehicle provided with 
a hydraulic power steering device, as described later. 

[0032] The road surface friction state estimating apparatus includes a 
steering torque sensor 1 1 which detects steering torque, a current sensor 12 which 

20 detects motor current, a steering angle sensor 13 which detects a steering angle, a 

vehicle speed sensor 14 which detects a vehicle speed, a lateral acceleration sensor 15 
which detects a lateral acceleration, a yaw rate sensor 16 which detects a yaw rate, 
and an electronic control unit (hereinafter, referred to as "ECU") 20 which estimates a 
road surface friction state using signals output from the sensors. 

25 [0033] The steering torque sensor 1 1 is disposed coaxially with a steering 

shaft. The steering torque sensor 1 1 outputs a sensor signal according to steering 
torque applied to a steering shaft, and supplies the sensor signal to the ECU 20. The 
current sensor 12 outputs a sensor signal according to the motor current of an electric 
motor which is used in the electric power steering device, and supplies the sensor 

30 signal to the ECU 20. 

[0034] The steering angle sensor 13 outputs a sensor signal according to a 
steering angle 9 P due to the driver's steering operation, and supplies the sensor signal 
to the ECU 20. The vehicle speed sensor 14 outputs a sensor signal according to a 
vehicle speed (a vehicle body speed) u, and supplies the sensor signal to the ECU 20. 
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[0035] The lateral acceleration sensor 15 outputs a sensor signal according 
to an acceleration in a lateral direction of the vehicle (a lateral acceleration), and 
supplies the sensor signal to the ECU 20. The yaw rate sensor 16 outputs a sensor 
signal according to a yaw rate which is a rotational angular speed around a center of 
5 gravity of the vehicle, and supplies the sensor signal to the ECU 20. 

[0036] FIG. 2 is a block diagram showing a functional configuration of the 
ECU 20. 

[0037] The ECU 20 includes a steering torque detecting portion 21 which 
detects steering torque; an assist torque detecting portion 22 which detects assist 

10 torque; a SAT estimating portion 23 which estimates SAT; a slip angle estimating 

portion 24 which estimates a slip angle; and a high pass filter 25 which performs high 
pass filter processing on the slip angle. 

[0038] The ECU 20 further includes a lateral force calculating portion 26 
which calculates a front wheel lateral force of the vehicle; a slip angle converting 

15 portion 27 which converts the front wheel lateral force into a slip angle; a low pass 
filter 28 which performs low pass filter processing on the slip angle obtained by 
conversion; a summing device 29 which sums two slip angles on each of which the 
filter processing has been performed; and a grip degree estimating portion 30 which 
estimates a grip degree; and a road surface |li value estimating portion 3 1 which 

20 estimates a road surface friction coefficient ja (hereinafter, referred to as "a road 
surface [i value"). 

[0039] The steering torque detecting portion 21 detects steering torque 
applied due to the driver's steering operation based on the sensor signal from the 
steering torque sensor 11, and supplies the detected steering torque to the SAT 
25 estimating portion 23. 

[0040] The assist torque detecting portion 22 detects the assist torque 
applied to the electric power steering device, based on the motor current based on the 
sensor signal from the current sensor 12 and a preset parameter (for example, a lead of 
a pinion, a lead of a ball screw, or an assist motor torque coefficient), and supplies the 
30 assist torque to the SAT estimating portion 23. The assist torque detecting portion 22 
may use a current command value output to a motor of the electric power steering 
device, instead of the motor current. 

[0041] The SAT estimating portion 23 calculates a sum of the steering 
torque detected by the steering torque detecting portion 21 and the assist torque 
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detected by the assist torque detecting portion 22, thereby estimating SAT generated 
between a road surface and a tire, while removing friction of a steering system. 

[0042] FIG. 3 is a chart showing an estimated SAT value with respect to 
the sum of the steering torque and the assist torque. A width between two straight 
5 lines shows a magnitude of a hysteresis characteristic due to the friction of the 
steering system. An inclination of each straight line is 1. 

[0043] FIG. 4 is a chart showing the estimated SAT value with respect to 
the sum of the steering torque and the assist torque, for describing a method of 
removing the hysteresis characteristic. 
10 [0044] When the sum of the steering torque and the assist torque is 0 and 

the slip angle is also 0 while the vehicle moves straight, hysteresis does not occur, and 
the estimated SAT value is 0. 

[0045] When steering is performed and SAT is generated, the estimated 
SAT value is calculated using an inclination K, with respect to the sum of the steering 
15 torque and the assist torque. More specifically, the SAT estimating portion 23 

performs a calculation according to the following equation (1) using a discrete logic. 
[0046] 

T SA1 (k + 1) = T SA7 (k) + K x ' (T DA (k + 1) - T DA (k)) -(I) 

[0047] In the equation, T SAT is the estimated SAT value, T DA is the sum of 
20 the steering torque and the assist torque. The inclination Kj is set to be smaller than 1 
in order to indicate that a variation in the estimated SAT value is small even when the 
sum of the steering torque and the assist torque varies due to Coullomb Friction or the 
like. 

[0048] When steering is performed, the estimated SAT value obtained 
25 according to the equation (1) reaches a point A in FIG. 4. Further, when the sum of 
the steering torque and the assist torque increases, the estimated SAT value increases 
according to a straight line which indicates a lower limit of the model, that is, the 
estimated SAT value increases according to an equation (2). 
[0049] 

30 T SA7 (k + 1) = T^k) + T DA (k + 1) - T DA (k) .-.(2) 

[0050] When steering is performed further, the estimated SAT value 
reaches a point B at which the steering angle stops increasing, and then the sum of the 
steering torque and the assist torque decreases, the estimated SAT value decreases 
according to the equation (1) using the inclination K,. In this region, the estimated 
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SAT value is set such that a variation in the estimated SAT value is small with respect 
to a variation in the sum of the steering torque and the assist torque. The estimated 
SAT value is thus set so as not to be influenced by the Coullomb Friction of the 
power steering device, or the like even if the steering force supplied by the driver 
5 varies to a certain degree when the steering angle is maintained while the vehicle is 
turning. 

[0051] When the estimated SAT value decreases from the point B to a 
point C, and then the sum of the steering torque and the assist torque increases, the 
estimated SAT value increases toward the point B according to the equation (1). 

10 Also, when the sum of the steering torque and the assist torque decreases from the 
point C due to a decrease in the steering angle, and the estimated SAT value reaches 
an upper limit of the model, the estimated SAT value decreases according to a straight 
line indicating the upper limit, that is, the estimated SAT value decreases according to 
the equation (2). Setting two types of inclinations in this way removes the hysteresis 

15 characteristic. 

[0052] FIG. 5A is a chart showing a time-dependent waveform of the sum 
of the steering torque and the assist torque during running on a road having a high |u 
value, and FIG. 5B is a chart showing a time-dependent waveform of the estimated 
SAT value during running on the road having a high |u value. The estimated SAT 

20 value in FIG. 5B is obtained by removing the hysteresis characteristic from the sum of 
the steering torque and the assist torque in FIG. 5 A according to the equations (1) and 
(2). By removing the hysteresis characteristic, a variation in the estimated SAT value 
when the steering angle is maintained, which is assumed to be caused by influence 
such as Coullomb Friction, can be almost completely compensated. 

25 [0053] FIG. 6A is a chart showing a relation of the sum of the steering 

torque and the assist torque with respect to a slip angle during running a road having a 
high ]i value and during running on a road having a low |a value. FIG. 6B is a chart 
showing a relation of the estimated SAT value with respect to the slip angle during 
running on the road having a high |i value and during running on the road having a 

30 low jlx value. In FIG. 6B, the estimated SAT value changes substantially linearly with 
respect to the slip angle both in the case of running on the road having a high \x value 
and in the case of running on the road having a low p value. Therefore, it is apparent 
from FIG. 6B that the hysteresis characteristic is removed. 
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[0054] The slip angle estimating portion 24 estimates a front wheel slip 
angle ot E [rad] which is a slip angle of a front tire, based on a steering angle 9 P [rad] 
and a vehicle speed u [m/s]. The steering angle 9 P is based on the sensor signal from 
the steering angle sensor 13, and the vehicle speed u is based on the sensor signal 
from the vehicle speed sensor 14. The front wheel slip angle ct E is represented by 
state equations (3) and (4), using dynamic characteristics of the vehicle movement. 

[0055] 
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[0056] In the equations, v is a lateral speed [m/s], r is a yaw rate [rad/s], u 
is a vehicle speed [m/s], c f is a front wheel cornering power [N/rad], c r is a rear wheel 
cornering power [N/rad], L f is a distance between a front axle and a center of gravity 
of the vehicle [m], L r is a distance between a rear axle and a center of gravity [m] of 
the vehicle, M is a vehicle mass [kg], I z is yaw inertia [kgm 2 ], and g h is a steering gear 
ratio. 

[0057] By making the equations (3) and (4) discrete using a sample time x 
so that they become functions of the vehicle u, the following equations (5) and (6) are 
obtained. 

[0058] 
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a E (k) = [1 L f ]x(k) I u(k) - d p (k)/g h - (6) 

[0059] In the equations, k is a sampling number. Also, A s and B s in the 
equation (5) are represented by the following equations (7). 
[0060] 
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[0061] The slip angle estimating portion 24 estimates the front wheel slip 
angle ct E by performing calculations according to the equations (5) to (7), and supplies 
the front wheel slip angle a E to the high pass filter 25 every sample time x. 

[0062] The high pass filter 25 performs high pass filter processing on the 
front wheel slip angle a E estimated by the slip angle estimating portion 24. When a 
steering neutral point moves during running on a bank road, the front wheel slip angle 
a E estimated by the slip angle estimating portion 24 contains a drift error in a low 
frequency region, but contains a signal component which does not have a phase lag 
with respect to the estimated SAT value in a high frequency region. Therefore, the 
high pass filter 25 removes the drift error in the low frequency region and extracts 
only the high frequency component which does not have a phase lag with respect to 
the estimated SAT value, by performing the high pass filter processing on the front 
wheel slip angle oc E . 

[0063] The high pass filter 25 is configured as a primary discrete filter. 
The primary high pass filter during a continuous time is represented by the equation 
(8) which is a transfer function. 

[0064] 

G H (*) = — — "(8) 
s + co b 

[0065] In the equation, co b is a breakpoint frequency. The high pass filter 
during a discrete time can be designed by transforming the equation (8), for example, 
by a Tustin transform. In the Tustin transform, when the sampling time is T and a 
time lead operator is z, a differential sign s is represented by the equation (9). 

[0066] 

s = ^zR ... (9 ) 

TXz + l) 

[0067] By substituting the equation (9) into the equation (8), the high pass 
filter during a discrete time is represented by the equation (10). 
[0068] 
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[0069] The high pass filter 25 performs the high pass filter processing on 
the front wheel slip angle ot E according to the equation (10), and supplies the 
summing device 29 with the front wheel slip angle a E on which the high pass filter 
processing has been performed. 

[0070] The lateral force calculating portion 26 calculates a front wheel 
lateral force F f which is a lateral force generated in a front tire, using a lateral 
acceleration g y and a yaw rate r. The lateral acceleration g y is based on the sensor 
signal from the lateral acceleration sensor 15, and the yaw rate r is based the sensor 
signal from the yaw rate sensor 16. The front wheel lateral force F f satisfies the 
equation of motion (11) concerning the lateral acceleration g y , and satisfies the 
equation of motion (12) concerning the yaw rate r. 

[0071] 

K^ + ™] =F/+Fr '" (n) 

I z ^= Lf F f -L r F r ...(12) 

[0072] In the equations, F r is a rear lateral force. The lateral acceleration 
g y is represented by the following equation (13). 
[0073] 

S,=£ + r» ...03) 

[0074] Based on the equations (11) and (12), the front wheel lateral force 
F f is represented by the following (14). 
[0075] 

F f = 2L ...(14) 

1 L f + L r 

[0076] The lateral force calculating portion 26 calculates the front wheel 
lateral force F f using the yaw rate r and the lateral acceleration g y according to the 
equation (14), and supplies the slip angle converting portion 27 with the front wheel 
lateral force F f . 
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[0077] The slip angle converting portion 27 converts the front wheel 
lateral force F f into a front wheel slip angle a T by dividing the front wheel lateral 
force F f supplied by the lateral force calculating portion 26 by a front wheel cornering 
power c f . More specifically, the slip angle converting portion 27 performs a 
5 calculation according to the following equation (15). 

[0078] 




10 [0079] The low pass filter 28 performs low pass filter processing on the 

front wheel slip angle ot x calculated by the slip angle converting portion 27. The 
front wheel slip angle a x calculated by the slip angle converting portion 27 contains, 
in the high frequency region, a variation component such as noise due to road surface 
disturbance and a phase lag. However, the front wheel slip angle a T contains the low 

15 frequency component which is not influenced by the road surface disturbance even 
when the vehicle runs on a bank road. Thus, the low pass filter 28 removes the 
variation component in the high frequency region and extracts only the low frequency 
component which is accurately calculated, by performing the low pass filter 
processing on the front wheel slip angle ot T . 

20 [0080] More particularly, the low pass filter 28 is configured as a primary 

discrete filter having the same breakpoint frequency as that of the high pass filter 25. 
The primary low pass filter during a continuous time is represented by the following 
equation (16) which is a transfer function. 
[0081] 

25 

G L (s) = -^S- -(16) 
* + co h 

[0082] The low pass filter during a discrete time is obtained by 
transforming the equation (16) by the Tustin transform, and is represented by the 
following equation (17). 
30 [0083] 

tohto±» ... (17) 

L (To) b + 2)z + Tco b -2 
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[0084] The low pass filter 28 performs the low pass filter processing on 
the front wheel slip angle cc T .according to the equation (17), and supplies the 
summing device 29 with the front wheel slip angle oc T on which the filter processing 
5 has been performed. 

[0085] The invention is not limited to a specific breakpoint frequency. 
However, the breakpoint frequency is preferably set so as to change in accordance 
with a road surface cant change speed in order to remove noise due to road surface 
disturbance. 

10 [0086] The summing device 29 calculates an integrated slip angle a x by 

summing the front wheel slip angle a E supplied by the high pass filter 25 and the front 
wheel slip angle a T supplied by the low pass filter 28. In other words, the summing 
device 29 performs a calculation according to the following equation (18). 
[0087] 

15 

a,(Z)= G„(Z)- a E (Z) + G L (Z)- a^Z) ». (18) 

[0088] The sum of the transfer function of the high pass filter 25 and the 
transfer function of the low pass filter 28 is 1. Therefore, when the same signal is 

20 input to the high pass filter and the low pass filter, and outputs of both the filters are 
summed, the original signal is restored. Accordingly, the summing device 29 can 
calculate the slip angle a, which is not influenced by a drift error or noise. 

[0089] FIG. 7 is a chart showing time-dependent waveforms of the 
estimated slip angle value a E , the converted slip angle value a T (i.e., the slip angle 

25 value obtained by conversion), and the integrated slip angle a, when the vehicle enters 
a bank road. In this experiment, the vehicle entered a bank road around 5 second, and 
then, the vehicle was steered to the left around 6 second. Owing to this vehicle 
movement, the estimated slip angle value a E became a negative value. When the 
vehicle runs on a bank road, the steering neutral point moves (in FIG. 7, the steering 

30 neutral point moves in a positive direction). Therefore, the slip angle used by the grip 
degree estimating portion 30 needs to be modified according to the shift in the 
steering neutral point. 

[0090] The estimated slip angle value a E estimated based on the steering 
angle value became 0 until steering was started even after the vehicle entered the bank 
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road (in an interval from 5 second to 6 second). The converted slip angle value <x T 
calculated based on the lateral force moved in the positive direction in accordance 
with the shift in the steering neutral point when the vehicle entered the bank road. 
However, the converted slip angle value a T contained much noise due to road surface 
5 disturbance. Meanwhile, it is apparent that the integrated slip angle ctj , which was 
obtained by integrating the signals using the filter processing, did not contain noise, 
and moved in the positive direction in accordance with the shift in the steering neutral 
point. 

[0091] The grip degree estimating portion 30 estimates a grip degree e 
10 which is equivalent to an excess of a frictional force in the lateral direction, based on 
the estimated SAT value T SAT estimated by the SAT estimating portion 23 and the 
integrated slip angle a! calculated by the summing device 29. 

[0092] The grip degree s is a ratio between the estimated SAT value T SAT 
and a SAT reference value T SATO , and is represented by the following equation (19). 
15 [0093] 

e = I™- -.(19) 
T 

SATO 

[0094] The SAT reference value T SATO is a SAT value of a linear model, 
20 that is, a model which is made linear using a slip angle 0, and is a SAT value in a high 
grip state, that is, a state where a grip degree is 1. When an inclination of the SAT 
value with respect to the slip angle at a zero point is Kq, the SAT reference value 
T SATO is represented as a product of the inclination K^, and the integrated slip angle a,. 
Therefore, the grip degree 8 is represented by the following equation (20). 
25 [0095] 

£ ^ISAT_ = T SAT ... (20) 

T S ato K 0 - a x 

[0096] The grip estimating portion 30 calculates the grip degree 8 using 
30 the estimated SAT value T SAX and the integrated slip angle <x„ according to the 

equation (20), and supplies the grip degree estimating portion 31 with the grip degree 



8. 
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[0097] The invention is not limited to a configuration in which the grip 
degree estimating portion 30 estimates the grip degree e using the aforementioned 
method. For example, the grip degree e may be represented using a function of the 
SAT reference value T SATO and the estimated SAT value T SAT . Alternatively, the grip 

5 degree s may be represented using a two-dimensional map of the SAT reference value 
T SATO and the estimated SAT value T SAT . 

[0098] FIG. 8 is a chart showing the estimated SAT value with respect to 
the slip angle during running on a dry asphalt bank road. A solid line indicates a case 
where the integrated slip angle otj is used as a slip angle, and a dotted line indicates a 

10 case where the estimated slip angle value a E is used as a slip angle. For reference, a 
dashed line indicates the SAT reference value T saxo . 

[0099] FIG. 9A and FIG. 9B indicate data measured in the experiment 
shown in FIG. 8. FIG. 9A is a chart showing a grip degree during running on a dry 
asphalt bank road, and FIG. 9B is a chart showing a time-dependent waveform of the 

15 steering angle. In FIG. 9 A, a solid line indicates a grip degree calculated using the 
integrated slip angle a l5 and a dotted line indicates a grip degree calculated using the 
estimated slip angle value a E . Further, when the calculated grip degree exceeds 1, the 
grip degree is set to 1 . 

[00100] Since the estimated slip angle value a E is influenced by the shift in 

20 the steering neutral point during running on a bank road, the estimated slip angle 
value ct E varied in the negative direction, as shown in FIG. 8. As a result, the 
absolute value of the estimated slip angle value cc E became small as compared with 
the absolute value of the SAT reference value T SAT0 . Also, the grip degree became 
small, as shown in FIG. 9 A. 

25 [00101] Meanwhile, the integrated slip angle a, is not influenced by the 

shift in the steering neutral point even during running on a bank road. In other words, 
as shown in FIG. 8, the relation of the estimated SAT value with respect to the 
integrated slip angle a Y substantially matched a straight line indicating the SAT 
reference value T SAT0 . Also, the grip degree became 1, as shown in FIG. 9A. 

30 [00102] According to the results of the experiment described above, in the 

case of a dry asphalt bank road, the grip degree constantly needs to be approximately 
1. However, the grip degree calculated using the estimated slip angle value oc E 
became a value smaller than 1, and it is apparent that this value is erroneous. 
Meanwhile, the grip degree calculated using the integrated slip angle a, constantly 



10 
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became 1. Therefore, it was confirmed that the grip degree was accurately estimated 
when the integrated slip angle a, was used. 

[00103] The road surface \i value estimating portion 31 estimates a road 
surface jx est value based on the grip degree e and the lateral acceleration g y when the 
grip degree s estimated by the grip degree estimating portion 30 becomes equal to or 
smaller than a predetermined reference value (for example, when the grip degree s 
becomes equal to or smaller than 0.5). The road surface |i est value is represented by 
the following equation (21). 

[00104] 



/'~=7T £ T---- (21) 



[00105] In the equation, g is a gravitational acceleration. Also, g fy is a front 
wheel position lateral acceleration, and is represented by the following equation (22). 
15 [00106] 

I z dr 

MJx_ ... (22) 

[00107] The estimation accuracy of the road surface ^ est value thus obtained 
20 is improved as the grip degree e decreases, that is, the grip degree 8 approaches a limit 
value. Accordingly, when the grip degree s becomes equal to or smaller than the 
predetermined reference value as described above, the road surface \x. estimating 
portion 31 estimates the road surface ^i est value according to the equations (21) and 
(22). 

25 [00108] FIG. 1 OA is a chart showing the grip degree during running on an 

ice road at a speed of 40 km/h, FIG. 1 0B is a chart showing the road surface |u est value 
(an estimated \x value), and FIG. 10C is a chart showing a time-dependent waveform 
of the steering angle. In FIG. 10B, a determination flag becomes 1 when the grip 
degree becomes equal to or smaller than the reference value, and becomes 0 when the 

30 grip degree exceeds the reference value. According to FIG. 10B, when the 

determination flag was 1, the estimated \x value constantly became approximately 0.3 
which is substantially equal to the road surface \i value of an ice road. 
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[00109] FIG. 1 1 A is a chart showing the grip degree during running on a 
pressed snow road at a speed of 40 km/h, FIG. 1 IB is a chart showing the road 
surface ji est value (the estimated value), and FIG. 1 1C is a chart showing a time- 
dependent waveform of the steering angle. According to FIG. 1 IB, when the 

5 determination flag was 1, the estimated ja value constantly became large as compared 
with the case of running on the ice road, and constantly became approximately 0.5 
which is substantially equal to the road surface ja value of a pressed snow road. 

[00110] FIG. 12A is a chart showing the grip degree during running on an 
asphalt road at a speed of 40 km/h, FIG. 12B is a chart showing the road surface |a est 

10 value (the estimated |u value), and FIG. 12C is a chart showing a time-dependent 

waveform of the steering angle. According to FIG. 12B, the grip degree constantly 
became equal to or larger than the reference value (0.5), and the determination flag 
constantly became 0. In other words, when running on a dry asphalt road at a speed 
of 40 km/h, there was a sufficient frictional force left between the tires and the road 

15 surface, though steering was performed in approximately the same manner as in the 
experiments shown in FIG. 10C and FIG. 1 1C. Thus, the grip degree s was higher 
than the reference value, and accordingly, the estimated |a value was not estimated. 

[00111] As described above, the road surface friction state estimating 
apparatus according to the embodiment can accurately estimate the grip degree s 

20 which is equivalent to an excess of a frictional force in the lateral direction, by 

calculating the SAT reference value T SAXO based on the integrated slip angle that has 
been extracted by the high pass filter 25 and the low pass filter 28, and calculating a 
ratio between the SAT reference value T SATO and the estimated SAT value T SAT 
estimated by the SAT estimating portion 23. In addition, the road surface friction 

25 state estimating apparatus according to the embodiment can accurately estimate the 
road surface ja value when the grip degree becomes equal to or smaller than the 
reference value. 

[00112] Particularly, the road surface friction state estimating apparatus can 
accurately estimate a road surface friction state such as the grip degree s and the road 
30 surface \i value, by making the estimated slip angle value a E pass through the high 
pass filter 25 so as to remove a drift error caused by the shift in the steering neutral 
point during running on a bank road. Further, the road surface friction state 
estimating apparatus can accurately estimate a road surface friction state without 
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being influenced by road surface disturbance, by making the converted slip angle 
value a T pass through the low pass filter 28. 

[00113] The breakpoint frequency of the high pass filter 25 and the low 
pass filter 28 may be a fixed frequency. However, the breakpoint frequency needs to 

5 be set to be equal to or higher than a road surface cant change speed, that is, a speed at 
which a cant of a road surface changes. The road surface cant change speed is 
proportional to a vehicle speed. Thus, the high pass filter 25 and low pass filter 28 
may be configured such that the breakpoint frequency becomes higher as the vehicle 
speed becomes higher. 

10 [00114] Since the road surface friction state estimating apparatus is 

provided with the high pass filter 25 and the low pass filter 28 thus configured, the 
apparatus can accurately estimate a road surface friction state even when the vehicle 
enters a bank road at a high speed. 

[00115] Also, the high pass filter 25 and the low pass filter 28 may be 

15 configured such that the breakpoint frequency becomes higher as a deviation between 
the estimated slip angle value cc E and the converted slip angle value a T becomes 
larger. This is because the deviation between the estimated slip angle value a E and 
the converted slip angle value a x becomes large in the case where the vehicle runs on 
a bank road, or in the case where the relation between the lateral force and the slip 

20 angle becomes non-linear. In such a case, it is preferable to use the converted slip 
angle value ot T which is based on the front wheel lateral force F f . that is not 
influenced by a change in the steering neutral point or the non-linear relation between 
the lateral force and the slip angle. 

[00116] Since the road surface friction state estimating apparatus is 

25 provided with the high pass filter 25 and the low pass filter 28, the apparatus can 

increase a proportion of the converted slip angle value a T based on the front wheel 
lateral force F f to the estimated slip angle value oc E based on the steering angle when 
using the converted slip angle value a T and the estimated slip angle value ot E to 
obtain the integrated slip angle ot I? according to the vehicle movement state. As a 

30 result, the road surface friction state estimating apparatus can accurately estimate a 

road surface friction state, for example, even when the vehicle suddenly enters a bank 
road, or the vehicle is suddenly brought into a spin state. 
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[00117] The invention is not limited to the aforementioned embodiments, 
and various modifications can be made to the design within the scope of the invention 
defined by claims. 

[00118] For example, in the aforementioned embodiments, the grip degree 
5 and the road surface |a value are estimated using the electric power steering device. 

However, a hydraulic power steering device may be used. When the hydraulic power 
steering device is used, it is possible to estimate the grip degree and the road surface |a 
value in the same manner as in the aforementioned embodiments, by measuring a 
hydraulic pressure and the like of the hydraulic power steering device, and detecting 
10 torque corresponding to the steering torque and torque corresponding to the assist 
torque. 

[00119] In the aforementioned embodiments, the high pass filter 25 and the 
low pass filter 28 are represented using the linear transfer functions. However, other 
functions may be used. 



